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L!x. \ Abstract: We investigate possible symmetry properties of the scalar sector of Multi- 

Higgs-Doublet Models, and, to some extent, the generalization of such models to gauge 

C^^ ' groups other than SU{2)l x U{1)y- In models where the C (charge conjugation) violating 

^^ ■ operator C is not present, the scalar potential is invariant under a group larger than the 

gauge group, 0(4) when the Higgs fields are doublets. If the Higgs fields develop aligned 
vacuum expectation values, this symmetry will break to an 0(3) subgroup, which in general 

k^ , is further broken by loop corrections involving the gauge bosons. Assuming such correc- 

H ' tions are small, the physical properties of the Higgs sector will approximately organize into 

representations of SO {3). If the vacuum expectation values of the Higgs fields are aligned in 
the direction of the C even fields, the mass spectra of the charged and C odd sectors will be 
degenerate. Moreover, if the Higgs fields develop a pair of non-aligned vacuum expectations 
values, so that the charge conjugation symmetry is spontaneously broken (but not the U{1) 
electromagnetic gauge invariance), a pair of light charged Higgs bosons should appear. 
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1. Introduction 

Future discoveries at the LHC may reveal a sector of scalar particles that is much richer than 
that of the Standard Model (SM). Some of the scalars may be responsible for generating 
the masses of fermions and the electroweak bosons [1], whereas others may be responsible 
for the dark matter [2, 3, 4, 5, 6, 7, 8]. 

It is natural to classify such scalars according to their properties under the SU{2) 
associated with the electroweak sector of the Standard Model. In order to be compatible 
with electroweak precision data, one usually considers only SU{2) doublets and singlets. 
Even these representations are severely constrained by the data [9, 10, 11, 12]. 

Going beyond one or two doublets [13, 14], one immediately has to face models having 
a large number of parameters. The structure of such potentials has been studied in [15]. 

Different doublets could be distinguished via their couplings to fermion fields. This 
idea is exploited in the so-called Model II version of the two-Higgs-doublet model (2HDM) 
[13], where one doublet couples to up-type fermions, and the other couples to down-type 
fermions. Another version of this idea is the one considered in Ref. [7], where each fermion 
or each family has its own Higgs field. 

We shall here consider instead the case when the different doublets can not be distin- 
guished (since we are not considering couplings to the fermions). Thus, we shall here study 



the symmetries of models with N doublets — it turns out that by setting a certain S'0(4)- 
violating parameter A''^-' to zero and assuming vacuum alignment, the spectrum simplifies 
considerably. In particular, a certain "custodial" 50(3) symmetry [16] leads to a degeneracy 
between the mass matrix of the C odd (or equivalently, CP odd) and the charged Higgs 
bosons. 

This possibility of a symmetry group of the scalar polynomial which is larger than 
required by gauge invariance was pointed out by Weinberg [17] many years ago. In the 
theories considered the extra symmetry was assumed to be a symmetry of any quartic (i.e. 
renormalizable) potential of the scalar sector. 

The Standard Model with its single Higgs doublet is an example of a theory were the 
most general scalar potential has an extra 0(4) symmetry not contained in the SU{2) x U{1) 
gauge symmetry. On the other hand, an extension of the SM with an extra Higgs doublet (as 
required if we want to introduce supersymmetry) or more, and with C symmetry, induces a 
most general scalar potential where the extra 0(4) symmetry is broken (before spontaneous 
symmetry breaking). In our notation, this extra symmetry is broken by the parameter (s) 
A'^', cf. eq. (2.8). Standard renormalizability instructs us to include into the Lagrangian, 
all terms allowed by the SU{2) x U{1) symmetry, hence it may be inconsistent to leave 
out terms proportional to A^^-* . The presence of such parameters will in general lead to an 
order A*-^-* tree-level breaking of the additional symmetry. We don't think this in principle 
is any different from having the symmetry broken by loop corrections. It becomes in any 
case a question of the magnitude of the perturbation. 

However, for renormalization it is only necessary to study the theory with unbroken 
symmetry, since the renormalization is not changed when f 7^ [26]. Hence, only quadratic 
counterterms have to be considered in the case of mass relations. Thus, in this case coun- 
terterms of the type A'^-* will be irrelevant. 

2. The NHDM potential and Lagrangian 

We define the A^-Higgs-doublet model, abbreviated NHDM, to be a system of N two-com- 
ponent complex scalar fields <I>i, <I>2, . . . , ^at, each with the same transformation property 
under SU(2)l x U{1)y as the Higgs field of the Standard Model, and with dynamics defined 
by the Lagrangian density 

C{X) = Y, [D^''^m{x)]HDf,'^ra{x)] " ^(^1, ^2, • • • , ^n), (2.1) 

m 

where V{^i,^2, ■ ■ ■ ,^n) is a potential that — in its most general form — is given by (2.5) 
below. The covariant derivative D^ is defined as 

D^ = a^ + igTiWl" + ig'YBf", (2.2) 

where Wj^ and B^^ are the SU{2)l and U{1)y gauge fields, respectively, and Tj = ^'^i 
are the generators of SU{2) with o"i,o"2,o'3 the Pauli matrices. Thus, our Higgs fields are 
labeled by two indices: The row index running from 1 to A^ is often written out explicitly as 
above, and an often hidden group index acted on by the gauge group. The latter are acted 



^m = 


^]n^m, 






-Dmn ^ 


U<^l^n + <^i^^) 


= Re(C<J>„) 


^Ba, 


^mn ^ 


-l{^l^„-^l^^) 


= Imi^l^n) 


^Ca. 



on by the matrices Tj in (2.2) (whose indices are also hidden). When written exphcitly we 
shah use Greek letters from the beginning of the alphabet. 

To write the most general gauge-invariant potential in a renormalizable NHDM in a 
compact way, we introduce a set of linearly independent^ hermitian operators invariant 
under local SU{2)l x U{1)y transformations (this is a generalization of the approach for 
the 2HDM in [18]): 



(2.3) 



Due to (anti-)symmetry under interchange of m and n we may impose the restriction that 
1 < 771 < n < iV, and introduce indices a, b, . . . labeling such pairs. An explicit invertible 
encoding is 

1 < a = a{m, n) = m + -{n - 2){n - 1) < -N{N -1)=M. (2.4) 

We let m{a),n{a) denote the inverse of this encoding. We will use the summation convention 
that repeated indices from the beginning of the alphabet are summed from 1 to M, and 
repeated indices from the middle of the alphabet are summed from 1 to A^. The most 
general potential in the NHDM thus becomes a linear combination of all different quadratic 
and quartic factors in the $„ (and ^m) which can be formed from Am, Ba and Ca'-^ 

yg($l, ...,'^n)= /i^^lm + f^^a^Ba + fl^^^Oa + A« 1„A„ + X^^b ^aBb 



+ >'ab CaCb + X^J^lAmBa + A(^) 1„C, + A^J BaCb, (2.5) 



where the "g" in Vg denotes "general". To avoid double counting we introduce the restriction 
m < n in the term involving Amn, and the restriction a < 6 in the terms involving X^fj , 
A^^ and A^^;, . We will not consider terms of degree higher than four, because these would 
destroy the renormalizability of the model [22]. From the hermiticity of the potential Vg 
all parameters fi and A in the expansion (2.5) must be real. Thus the number of free real 
parameters in (2.5) is 

Ntot = N + 2Ar + ^N{N + I) + AfiAT + 1) + 2NAr + Af^ = ^N\N^ + 3), (2.6) 

which for A^ = 1 gives us the 2 parameters of the Standard Model (//^ and A). N = 2 gives 
us the usual 14 parameters for the 2HDM. There are 54 parameters for A^ = 3 and 152 
parameters for A^ = 4. 

This counting does not take into account the fact that we may make SU{N) row 
transformations on the fields <!>„ to eliminate some terms in (2.5). One possible choice is 



^There are no linear relations between the operators in (2.3). However, they are algebraically dependent 
when N > 3, being restricted by (A'' — 2)^ polynomial equations of 8'th order in the fields. 

^Since (^la^^A ($J„o-^$„) = - (^l^i) i^L^n) + 2 (^l^„) {^l^^^i), other quartic invariants may 
be expressed by those chosen. 



(2) (3) 

to transform the quadratic terms into a diagonal form, i.e. so that fia = fJ-a = 0. This in 

general leaves a matrix of A^ — 1 independent diagonal phase transformations (such that the 

(5) 
determinant is unity). We may for instance use it to transform all A^^ with n{a) = m{a) + 1 

to zero. This reduces the number of parameters by iV^ — 1, i.e. to N^^^ = 2 {N^ + N"^ + 2) , 

yielding 11 parameters for A^ = 2 (in agreement with Barroso et. al. [14]), 46 parameters 

for A^ = 3, and 137 parameters for A^ = 4. 

2.1 The most general C-invariant NHDM-potential 

The charge conjugation operator C is a linear operator which leaves complex constants 
unaltered, but maps fields onto their hermitian conjugate; C{z<^m) = z^^ln, where z is a 
complex number.^ Then C{Ca) = —Ca, in contrast to C{Am) = Am, and C{Ba) = Ba- 
We obtain a C-invariant potential by leaving out all terms which are odd in Ca, i.e., terms 
involving fi^a\ \ma, and A^J. There are M + NM + M"^ = \N{N - 1){N'^ + N + 2) such 
terms, leaving 

Nc = -N{N^ + 5N + 2) (2.7) 

free parameters for the general renormalizable C-invariant potential, 

+ X^^daCh + X^JllAmBa. (2.8) 

For A^ = 1 we get the usual 2 parameters of the standard model, the Higgs potential of which 
is automatically C-invariant. For A^ = 2 we get the usual (see e.g., [18]) 10 parameters. For 
A^ = 3 we get 33 parameters, and for A^ = 4 we get 86 parameters. 

This counting does not take into account that we may make 0{N) transformations on 

the row of fields <&„ to eliminate some terms in (2.8). A natural choice is to transform 

(2) 
the quadratic terms to diagonal form, i.e. so that Ha =0. This reduces the number of 

parameters by 2-^(-^ ~ 1)) i-^- to A'^ = jN{N^ + 3N + 4). This gives 9 parameters for 
A'' = 2 (in agreement with [14]), 30 parameters for A^ = 3, and 80 parameters for A^ = 4. 

The difference A'phases = A't'ot - ^c = jN^i^'^ - 1) -^ 1 - A^ counts the number of 
genuine C-violating parameters in Vg (in agreement with Branco et al. [19]). 

2.2 Symmetries of A, B, C and C^ 

For generality we here consider k (rather than 2)-component fields, i.e. with SU{k) x U{1) 
as gauge group. To make it simpler to explore all possible symmetries we express the field 
$m in terms of its independent real (hermitian) components, <I>m = ^m + i0m- Define 
2k X 2k matrices 



^This definition assumes that we for some reasons have decided on a decomposition of all fields into their 
real and imaginary parts. It is not invariant under complex transformations of the fields, see e.g. [19]. 



where the subscript k mdicates the hnear dimension of the submatrix involved. The complex 
scalar product between two fields $m and <^„, invariant under unitary {U{k)) transforma- 
tions, can be expressed in terms of two real bilinear forms^ 

Re{^l<^n) = B„,n = (C 0^)^(0")= <^n + ©m©"' (2-lOa) 

Imi^l^n) = Cmn = «, 0™) ^ ( J = ^™0« - 0-^- (2.10b) 

The first is the Euclidean dot product between 2/c-component real vectors, the second is 
the Poisson bracket (symplectic form) of the same quantities viewed as coordinates and 
momenta of 2/i;-dimensional phase space. The quantities in (2.10) are individually invariant 
under transformation groups larger than U{k). The first form B (with A as a special case) 
is invariant under the 0{2k) group of real orthogonal transformations, 

(:;)-o(*»). oro-T, ,.n) 

the second form C is invariant under the Sp{k, R) group of real symplectic transformations. 

In this formulation the charge conjugation operator C discussed above can be represented 
as a particular 0{2k) transformation when acting on the fields ^„ and 0„ 



C = '• 7 . (2.13) 




Considering infinitesimal transformations, O = T + eL + C(e^), S =X + eM + 0{e'^), the 
conditions (2.11) and (2.12) become 

L^X + XL = 02fe, M^J + JM = 02k- (2.14) 

Thus L must be a 2A; x 2k antisymmetric real matrix; there is a set (Lie algebra) so{2k) of 
2k'^ — k linearly independent such matrices. Writing out the condition for M in terms of 
k X k submatrices we find that it must have the form 



M=\^ .rh B = B^,C = C^. (2.15) 




There is a set sp{k) of k'^ + k{k + 1) = 2k'^ + k linearly independent such matrices. The 
infinitesimal transformations of the original U{k) are the intersection of the sets so{2k) and 
sp[k). I.e., the matrices of the form ( „ . ), with A = —A and B^ = B. There are 



-B A 



^k{k — 1) + ^k{k + 1) = /c^ such linearly independent matrices. 



^{^m,^n) = $^ • $„ = Bmn + iC„ 



The symmetries of C^ The form C^ (or more precisely CmnCm'n') has a bigger sym- 
metry group than the form C. Still, we will see that such operators (forms) will violate the 
full 0(4) symmetry we can assign the rest of the Lagrangian. In analogy with (2.12), C"^ 
symmetries are given by 

which can be collected in a set 

P{k, M) = {5 G GL^kms'^JS = ±J}, (2.17) 

which we in Appendix A show is a Lie group. 
The component 

p-(fe,M) = {S& GL^kimS^JS = -J}, (2.18) 

consists of matrices with determinant 

det(p-(A;,IR)) = (-1)'=, (2.19) 

as shown in appendix A. The group P(/c,IR) = S'p(/c,M) U P^(A;,IR) will have the same Lie 
algebra as Sp{k,'M.), since the new component P~(/c,IR) is not connected with the identity. 
This is manifested by the equation corresponding to eq. (2.14), 

J + e{M^J + JM) = ±J. (2.20) 

not having any solution for the case of a —J on the right side, see appendix A for a proof. 

The most general 0(2/c)-symmetric potential We can conclude that the most general 
0(2/c)-invariant potential can be written 

Vb(2fc)(^l, . . . , ^7V) = ^^^^m + k^'-a^Ba + A^ l^^^ 

+ X^^^Bah + Xti^mBa, (2.21) 

since we have seen that operators not containing any factor C are 0(2/i;)-invariant. We 
obtain a 0(2A;)-invariant potential by leaving out terms proportional to A^^ from the C 
invariant potential Vc (2.8). The number of terms in Vo{2k) is then No{2k) = ^C ~ -A^^ = 
\N{N^ + 5iV + 2) - \N^{N - 1)2, giving 

Nom = \N{N + lf (2.22) 

free parameters for the general renormalizable 0(2/c)-invariant potential Vo{2k)- Foi" A^ = 1 
we get the usual 2 parameters of the standard model, the Higgs potential being automati- 
cally 0(2A;)-invariant. For A^ = 2 we get the usual 9 parameters, one parameter less than 
for the C-invariant potential. For A^ = 3 we get 24 parameters, and for A^ = 4 we get 50 
parameters. 

This counting does again not take into account that we may make 0{N) transformations 
on the row of fields <!>„ to eliminate some terms in (2.21). We may transform the quadratic 
terms to diagonal form, so that fXa =0. This reduces the number of parameters by 
\N{N - 1), i.e. to A^o(2fc) = \N{N'^ + N + 2). This gives 8 parameters for A^ = 2, 21 
parameters for A^ = 3, and 44 parameters for A^ = 4. 



2.3 Symmetries of the NHDM potential 

Since the NHDM-potential Vg is constructed from the invariants (2.3) the symmetries of the 
latter are reflected in the symmetries of the former, but in a manner depending on details 
of the construction: 

1. If Vg depends only on the Am's, i.e. if only the parameters fim and Amn are nonzero, 
then the symmetry group of Vg is at least^ ^rn=i 0{2k), since we can make indepen- 
dent transformations on each $m- 

2. If Vg depends only on the Am^s and the B^s, i.e. for a C-invariant theory (2.8), where 

(s) 
in addition the parameters A)jjj = 0, then the symmetry group of Vg is at least 0{2k). 

(2) (2) 

It may contain several such factors if some of the parameters Ha and A^^^ vanish. To 

(2) 

analyze this we partition the $rn's into sets: If a parameter Ha is nonzero, then the 
fields ^m(a) and ^n(a) belong to the same set, with m{a) and n{a) denoting that m 

(2) 

and n are contained in a. If a parameter A^^ is nonzero, then the fields ^m(a) and 
^n(a) belong to the same set, and the fields ^m(b) and ^n(b) belong to the same set. 
With this partitioning into a maximal number of sets we may make one independent 
0{2k) transformation for each set. 

^ (3) (3) 

3. If Vg depends only on the C^'s, i.e. with only the parameters fia and X^^ being 

nonzero, then the symmetry group of Vg is at least Sp{k, M). If we (in the same manner 
as above) can partition the fields into several sets, then we may make independent 
Sp{k, M) transformations on fields belonging to different sets. However, since the 
additional symmetries in this case fail to be symmetries of even the zero'th order 
kinetic term (2.30), their significance is uncertain. 

4. With all parameters arbitrary the symmetry group of Vg is just the original SU{k) x 
C/(l) gauge symmetry. 

In this work we will pay special attention to the second scenario, with k = 2. 

2.4 Symmetries of the kinetic terms 

We now turn to the (global) symmetries of the kinetic terms of the Lagrangian, 

N 

K = Y,[{d^ + G>')^n{x)]^[{d^, + G^)^n{x)], (2.23) 

n=l 

with 

G^" = igTiW^ + ig'YB^. (2.24) 

Let Ki denote the terms of the i'th order in the gauge fields. 



^It could possibly be larger, since there might be additional row symmetries transforming fields $„ with 
different m into each other; such symmetries would require special relations among the parameters /im and 



A 



(1) 



Consider the transformation of the kinetic terms linear in the gauge fields, Ki, under 
the map p defined in appendix B. We can then write^ 



N 



n=l 

N 

n=l 

TV 

= Y^ 9^Cr^$„ + <^l{-rnd,^n (2.25) 

n=l 

where the subscript k in ($n)fc indicates this is the usual complex Higgs k-p\et (in the case 
k = 2 the usual complex Higgs doublet), while $„ is the 2k dimensional real vector (2.26), 

*" = (e:)- P'^^' 

where {^n)k = ^n + ^©m and where we also use eqs. (B.2), (B.3) and (B.6). 

In eq. (2.25) we have applied the transformation p on the gauge terms G^ defined in 

eq. (2.24), 

p (G^) = T^", p (G^^) = -T'^, (2.27) 

where T then reads 

' -[gWi^ + g'YB^ gW^ J' ^'''"^ 

with W^ = ^- W^ Tf, summed over the set of real symmetric generators Tf of SU{k), and 
W^ = i^^iWJ^T^', summed over the set of imaginary antisymmetric generators T". For 
k = 2 the two sets are respectively {2^^) 2*^^} ^^'^ {2'^^}' 

Third, we consider the kinetic terms quadratic in the gauge fields, 

N 



K2 = Y,i^n)iG^^G^{^n)k 
ra=l 

N 

= Y.Pii^n)l)p{G^^)p{GM{^n)k) 
n=l 

N 

= -Y.'^nT^'^n. (2.29) 



n=l 



^Disregarding so-called Schwinger terms — here terms proportional to i[dn(f>{x) , (fi{x)] for a scalar field 
—or, alternatively, reasoning classically. 



The symmetries of Kq When we first ignore couplings to the gauge fields the remaining 
terms can be written 



N 



Ko = Y.Y. (^''*«- ^M*n« + 9^®na df^Ona) , (2.30) 

n=l a=l 

where the group index a labels the k components of $„ = ^„ + iG„. This term is invariant 
under rotation of all components {^na,Qna} into each other. I.e., the symmetry group of 
Kq is 0{2kN). The connected part of this group is S0{2kN), whose generators are all real 
antisymmetric matrices, Lmn,ai3 = —Lnrn,i3a (i.e. L = —L, where transposition refers to 
both sets of indices). 

The symmetries of Kq and Ki: Next, consider the terms linear in the gauge fields 
again, cf. eq. (2.25), 

^i = E((^MK,0^))r('Q;;V(^Le^)ra,(^^^')V (2.31) 

with T given in eq. (2.28). 

Consider now an infinitesimal transformation 6^rn,a = Lmn,ai3^n,i3, and T denoting 
the 2k X 2k antisymmetric matrix in equation (2.28) (in group indices a, /3 — in addition it 
is proportional to the N x N unit matrix in row indices). The requirement that this is an 
infinitesimal symmetry transformation for Ki is that L T+T L = 0. Or, when we restrict 
L to be antisymmetric so that it also is an infinitesimal symmetry transformation for Kq, 

J^mn,al3 '13^ ~ 'afS ^mn,l3^ = 0. (2.32) 

In order to determine the allowed structure of L, we expand these matrices into terms of 
definite symmetries {L^^' symmetric, and L^'^' antisymmetric) in the mn indices: 

^mn,al3 ^ / ^ ['^afi^mn "T -^afi J-' mn ) i (Z.oo) 

with S (symmetric) and A (antisymmetric) restricted by the constraint (2.32). The sum 
runs over all possible combinations of allowed matrices . We next note that the antisym- 
metric matrices T can be expanded in the set 



r={Tti,T:j,j} = ir: z]A L: \A " . • (2-34) 




By substituting (2.33) into (2.32) we are led to search for the set of 2k x 2k real matrices S 
and A which commute with T for arbitrary values of the fields W-^ and B^^. It is sufficient 
to verify that this property holds for all elements of the set T- Let 

^ = ( v" v'' ) ' ^ = -5' °^ ^ = ^- (2.35) 

\ A21 A22 / 



■^Without the restriction (2.32) there would be \k{2k + l)N{N -l) + \k{2k-l)N{N + 1) = kN{2kN-l) 
independent terms, equal to the number of generators of SO{2kN). 



Requiring commutativity [see eqs. (2.32) and (2.33)] with the three types of matrices in T 
we obtain the conditions 

















(2.36a) 


iiij = Jj A22, 


X22Tf 


= T!Xu, 


X12T! = 


-Ti X21, 


X21T! = 


-T!X,2, 


(2.36b) 



Xii=X22, Xi2 = -X21. (2.36c) 

Using (2.36c) we find that Xn and X12 must commute with ah matrices T^, Tf (assumed 
to form an irreducible representation). By Schur's lemma they must then be proportional 
to the k X k unit matrix, so that S (xX and A (x J . Thus, the Lie algebra of the symmetry 
group of Kq and Ki consists of elements of the form 

L„,n=TL^± + JL^±. (2.37) 

This is the Lie algebra of U{N) written in real variables. 

The symmetries of Kq and Ki in the hmit 51' — )■ 0: A more interesting situation 
arises if we remove J from the set T, as would apply to the limit g' — )■ 0. Then we still find 
that Xii + X22 and X12 — X21 must commute with all matrices T-^, Tf, and hence must be 
proportional to the unit matrix. Further, Xn — X22 and X12 + X21 must commute with 
all Tf", but anticommute with all T?. For k = 2, i.e. for the gauge group SU{2)l x U{1)y 
in the limit g' — t- 0, we find that nonzero solutions of (2.32), 

Xn - X22 a e = i ^2, X12 + X21 (X e, (2.38) 

are possible [see eqs. (2.33) and (2.35)]. This means that the possible antisymmetric ma- 
trices A may be any linear combination of matrices from the set 

J \ , (2.39) 

where the 2x2 matrix e was defined in eq. (2.38). The set Q is a, basis of generators for 
SU{2). Thus, the Lie algebra of symmetry generators for Kq and Ki in this case consists 
of elements of the form 

L^n=TL(:,l+Y,AL(;^l (2.40) 

Aeg 

allowing all possible symmetric N x N matrix L^^' for each A. There are ^N{N — 1) + 
^N{N + 1) = 2N'^ + N independent terms, equal to the number of generators of the N x N 
quaternionic symplectic group Sp{N). The generators (2.40) generate Sp{N), where the 
elements of Q act as the quaternions i, j and k. 

The results above were again found under the assumption that the fields Wj^ are ar- 
bitrary, and kept constant under the transformation. Combined SU{k) transformations of 
the Wj^ and the ^m fields still remain a symmetry. For k = 2 this symmetry is enhanced 
to (at least) Sp{N) x S0{4) as g' -^ 0.^ (In the case g' ^ it is U{N) x SUik).) The 
generators for 5*0(4) are the 3 generators in Q plus the 3 generators for the SU{2)l gauge 
group (written in real form). 

*The SO(4) symmetry cannot be extended to a 0(4) symmetry [27]. 
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The symmetries of K2 Next, consider the terms quadratic in the gauge fields cf. eq. (2.29), 

i^2 = -EK'0n)^'(el)- (2.41) 

As in the symmetry analysis of Ki we want to find all matrices X such that XT^ = T^X. 
All matrices X which commute with T will fulfill this criterion (since T^ can be expanded 
in a set which consists of products of all possible pairs of matrices from T). Therefore, the 
symmetries of Ki are also symmetries of i^2- 

3. Spontaneous symmetry breakdown 

In this section we return to the case of /c = 2, i.e. with SU{2)i x U{1)y as the gauge group 
and a row of N Higgs doublets <!>„. Note however that many of our considerations have 
straightforward generalizations to k > 2. 

As for the Standard Model, the potential Vg of equation (2.5), or Vq of equation (2.8), 
may acquire its minimum at nonzero values of the scalar fields, (<l>)o = ^^^\ where <I> 
(without a lower index) refers to the whole set of fields <!>„. This point, (^^^', will belong 
to one or more manifolds of equivalent minima related by the symmetries of the potential. 
One may use these symmetries to transform ^^^' to a particular form. A possible one is to 
require for (^\ that only its lowest real component is nonzero. This can always be achieved 
by an SU{2)l x U{1)y gauge transformation. Next, the upper component of $2 ^^"^ ^^ 
made real by the remaining U{1) gauge transformation which keeps <I>^ unchanged. Then 
one has no gauge freedom left to change <!>„ ' for n > 3. However, it was shown by Barroso 
et al. [23] that a sequence of unitary row transformations can shift the vacuum expectation 
values to the first two fields of the row only^ , for instance (when written in complex form) 

<' = („!)• *™ = („;:.). *!.»'= tor „> 3, (3.1) 

with vi, U2, V2 and 5 real. The special case U2 = (J = is usually referred to as vacuum 
alignment, in which case we may also transform V2 to zero by an orthogonal row trans- 
formation involving only $1 and ^2- This is known as the Higgs basis [24]. However, for 
other purposes it may be more convenient to adopt a "democratic" basis in which the lower 
component of all (or most) fields ^^ have a nonzero real expectation value. It is related 
to the Higgs basis by an orthogonal row transformation which preserves the form of C and 
C/(l) electromagnetic gauge transformations (the latter preserving the definition of electric 
charge) . 



®One may collect the quantities $ma (q = 1,2) into two A'^-coniponent complex vectors <i>'^' and $"'. 
By a U{N) row transformation one may first rotate <brn so that only the component <t[ is nonzero, with 
<i>j real. There is a group of U{N — 1) transformations preserving this condition; this may be used to 
transform $}„ so that only the components 9\ , $2 ^■'^ nonzero, with $2 real. One cannot do better 
due to the existence of four real U{N) invariant parameters in H^'^'H, ll"!'^'!!, and $(-'^)t i'^'. But there 
remains a U{N — 2) group of transformations preserving this condition which can be used for other purposes. 
For a SU{k) x f/(l) gauge group one may generalize this procedure to k vectors $'■-'', j = 1 . . . k. 
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Assume now the case of vacuum alignment and a potential Vq^^-j which is 0(4) invari- 
ant. Then the existence of the vacuum expectation values ^^^' will break the (explicit) 
symmetry down to 0(3), with the consequence that the Higgs boson particle spectra and 
other physical properties will organize themselves into multiplets of 0(3) (broken by per- 
turbative corrections in g'). The number of broken symmetry generators is 3 whether we 
consider the symmetry broken from U{2) to U{1) or from 0(4) to 0(3); this leads to the 
existence of 3 Higgs ghosts and no extra Goldstone bosons^'^. 

3.1 Mass-squared matrices 

To make these statements slightly more explicit, as needed for calculation of the zero'th 
order (in g and g') particle masses, we expand the potential around ^^^' to second order. 
There are no first order terms since we are expanding around a minimum. The matrix of 
second derivatives is the mass-squared matrix M^^^o . It is restricted by symmetries to have 



a block diagonal form in the group indices a, (3. We use coordinates where <I>^ 
is expressed in terms of four real fields, 



$. 



vl> + iB 



9ml + Wm2 
Vm ~r f]m i l<Pm3 



$(0) 




^m+ie. 



(3.2) 



It is now convenient to represent these on real form as 

/ 4'ml \ 




Vm ~r ?7»i 
0m2 



(3.3) 



\ (Pm3 J 



We have the expansion 

y($(0) + /^$) 



{V)o + 



d^v 



d^mp d<^r. 



A^rr,pA<^na + 0{A<^^ 



(3.4) 



where $mp denotes one of the four possibilities 



^ml ; 'Hm ) 9m2 j 



!>m3 



, and the subscript 



indicates that a quantity is evaluated at $ = ^^^' . Now a set of generators for 5*0(4)^^ is 

/O o\ 



Ji 



J2 



■h 



/O o\ 

0-10 

1 
\0 0/ 

/O -l\ 




\l o/ 

^°This remains true for general values of fc > 2; a set of aligned vacuum expectation values will break 
U{k) to U{k — 1) and 0{2k) to 0{2k — 1). The number of broken generators is 2fc — 1 in both cases. The 
situation is different if we have two broken real directions, as in equation (3.1) with U2 = but 5^0. Cf. 
section 3.3. 

^^Equivalently 0{4): SO{4) and 0(4) have the same Lie algebra and hence the same generators. 



•h 



I 1 o\ 
-10 

\ 0/ 

/O o\ 

1 

\o -1 o/ 



J5 



Jfi 



0-1 



\0 1 0/ 

/O -1 o\ 

10 

\0 0/ 



(3.5) 
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where Ji, J2, J3 will transform the vacuum expectation value ^^^' , while J4, J5, Jg leave it 
unchanged, cf. Eq. (3.3). In terms of these, the broken generators of the SU{2) x U{\)y 
gauge group, written in real form by the transformation p defined in (B.l), are 

\<y^^\{J2 + -h). V^i(Ji + J4), A(l-a3)^J3, (3.6) 

and the unbroken U{1) (electromagnetic gauge) generator is 

^(l + a3)^J6. (3.7) 

If now V is invariant under 

A(/>ml -> -A(/)mi, Ar/m -> A??™,, A(/)m2 -> -A(/>m2, A(/)m3 -> A0m,3, 

there can be no terms in (3.4) mixing the sets {/S.(j)i, /S.(j)2} and {At]^, A(/>m3}. If V in 
addition is invariant under 

A(/)„i -^ A(/)„2, Ar/„ -> A??™, A(/)„2 -> -A0m,i, A(/)„3 -> A(/)„3, 

there can be no terms in (3.4) mixing A0mi and /S.(j)m2, and we must have 

d'^V \ / d'^V \ _ , . ,^ ^, 

= ^ch,mn- (3-8) 



</^ml (^</>nl / \ ^9m2 UVn2 / q 

We refer to this as the charged mass-squared matrix. The transformations considered 
generate a Z4 subgroup of the U{1) gauge group generated by Jg, assumed to be a symmetry 
of V . We have formulated it this way as a reminder that invariance under discrete subgroups 
may be sufficient to impose useful restrictions on the mass matrices. 
If V is invariant under C transformations, 

Ar/„ -> A??m, A(/)m3 -> -A0m3, (3.9) 

(irrespective of how we define C to operate on the charged sector, f.i. A^^i — > A(/)mi, 
A0m2 — ^ — A(/)m2) there can be no terms in (3.4) mixing Ar/m and A(j)m3- Thus the neutral 
mass-squared matrix decomposes into two more blocks, a C even and a C odd one. 

If y in addition is invariant under the transformations 

A0ml -^ A0m,i, Ar/m -^ A??™,, A(^m2 -> A0m,3, A(^m3 -^ A0m,2, (3-11) 

which generate a Z4 subgroup of the S0{2) symmetry group generated by J4, we obtain 
the relation 

^C-,mn = ^ch,mn- (3.12) 

This explicitly demonstrates mass degeneracy between the charged and the C odd sec- 
tors [25]. Especially, if the potential is 0(4)-invariant (2.21), that is, we have a C-invariant 
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theory where also the parameters A^^ = [the latter implies (3.11)]^^, the above symme- 
try criteria for mass degeneracy are valid. Moreover, since the the renormalization is not 
changed when the Higgs fields aquire a vacuum expectation value [26], we won't get any mass 

renormalization counterterms from the quartic operators. So even though 0(4)-violating 

(3) 
quartic terms proportional to A^j, cannot alone be prohibited by any discrete symmetry 

imposed on the NHDM Lagrangian [27], they won't be showing up as counterterms when 

renormalizing the masses. Hence, the mass degeneration (3.12) will only be broken by loop 

corrections involving gauge bosons, since we get an exact S0{3) symmetry when g' = 0. 

With g' ^ and hence with an approximate 50(3) symmetry, the mass differences of the 

charged and CP-odd sectors will be of order 0{g"^) oc 0{e^). 

On the other hand, the 50(3) symmetry between OP-odd and charged sectors could 

also be broken by counterterms of the type X^^C"^, even though these terms are set to zero 

in the original potential, if we are considering scattering processes and not mass relations. 

3.2 The Higgs ghosts 

Let A<I> be chosen so that (^^^'+e/S.(^ + 0{e'^) is a family of minima related by the symmetry 
of the potential V ^ 

^ -y($(0) + eA$) = 0, (3.13) 



to first order in e. By differentiating this relation with respect to e and then setting e = 
we find 

which reflects the fact that the matrix M^^^o has zero eigenvalues with corresponding 
eigenvectors A$„^. We may take the latter to be A$'*' oc Ji^^^' for i = 1,2, 3. Normalized, 

A$« = iv„„ 0, 0, 0)^ /a, A$(^) = (0, 0, v„„ of /a, A^^) = (0, 0, 0, v„,f /a, (3.15) 

with a^ = ^„Wm- The massless excitations in these directions correspond to a triplet of 
Higgs ghosts. There will be A^ — 1 additional 50(3) triplets of excitations in directions 
orthogonal to the ghosts. They correspond to physical particles. There will also be A^ 
50(3) singlets, transforming evenly under C, corresponding to physical particles. In the 
case of A'' = 2, the triplet is {H~^ , H~ , A), whereas the singlets are h and H [13]. 

3.3 Non-aligned vacuum expectation values 

We have assumed vacuum alignment in much of the previous discussions of this section. 
The phenomenologically most realistic deviation from this case is that we have a situation 
with two (real) broken directions, as in (3.1) with ii2 = but (5 7^ 0. This corresponds to a 
situation which preserves the U{1) electromagnetic gauge symmetry, and its corresponding 
definition of electric charge, but where the C symmetry is spontaneously broken. In this 
situation the Mj^ ^^ mass-squared matrix remains in block form, but the C even and odd 
excitations may mix to give a 2A^ x 2A^ mass matrix for the neutral particles. One of the 



^For supersymmetric theories we typically have X^^ 7^ 0. 
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excitations will be massless, corresponding to a neutral Higgs ghost. We assume again that 
V is invariant under 0(4) transformations. The explicit symmetry is now broken down to 
0(2) ~ U{1), so that 5 generators are broken. As before, 3 of these will generate excitations 
which correspond to the Higgs ghosts; the remaining 2 will correspond to nearly massless 
charged pseudo-Goldstone^'^ bosons (massless to zero'th order in g'). 

To analyze the situation we again write $ = $(0) -)- $' in terms of real fields, 



<^m = ^^^ + {4>mi,rim,4>m2,Xmf with <^^^^ = (0,i;„,0,m 



T 
mj 



J4 and J5 are now also broken by the vacuum expectation values. Acting with the broken 
generators on ^^^' one finds five eigenvectors of the mass matrix with zero eigenvalues, 
A<I>(*) oc Jj$(°). After normalization 

A^g^ = {v^,0,0,0f /a, 
A<^g^ = {0,0,-vm,0f/a, 



A$(^) = (0, Wm, 0, -Vmf /vV + 62, (3.16) 

A$W = (o,o,«;„, 0)^/6, 
A$(^) = («;„, 0,0, 0)^/6, 

where a^ = X^^i'm ^^'^ ^^ — Ylm'^^rn- These eigenvectors are normalized, but they are not 
necessarily orthogonal to each other. Their nonvanishing inner products are 

A$W,A$(5)) =-('a$(2)^A$W) = ^Y^v^Wm^ cos l}. 

m 

Here | sini?| > 0, since the vacuum expectation values by assumption are non-aligned. Thus, 
the orthonormalized eigenvectors corresponding to the Higgs ghosts can be written 

^™^ = r^^na i^n.,0,-v^,0f = — ^= Acl>(^) + ^=L= A<^(^), (3.17) 

va"^ + V a + V a"^ + o 

where Hm oc Gi^^^\ Gi denoting the SU{2) generators as given by the map (3.6). The 
two eigenvectors corresponding to the Goldstone modes are orthogonal to those above, 

G^^) = ^ \a (a-^^^^ + cos T? A$(2)\ _ ^ (^^(i) _ cos ^ A$(^) 

sin T? V a^ + 5^ L V / V 

G(2) = ^'^ f-a (a^^^^ - cos 1? A^^i)) + b ("a^^^) + cos i? A$(^) 

sin T? V a^ + 6^ L V / V 



(3.1J 



^^Pseudo-Goldstone bosons stems from broken generators of the extra 0(4) symmetry of the potential, 
while Higgs ghosts per definition is generated by the broken generators of the gauge symmetry (which of 
course is a symmetry of the whole lagrangian). The pseudo-Goldstone bosons acquire small masses from 
radiative corrections, and are hence not massless in all orders of perturbation theory, like true Goldstone 
bosons. True Goldstone bosons are, in contrast, generated by the spontaneous breaking of a symmetry of 
a total lagrangian, not only a potential. 
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They have been orthonormahzed. We note that the normahzation constant becomes infinite 
in the hmit of ahgned vacuum expectation values, siwd — )■ 0. We recah that the set 
\^H''^' jH^"^' , H^^' ,G^^' ,G^'^'\ are just numerical eigenvectors of the mass-squared matrix. 
The corresponding zero mode fields are the quantum fields obtained by projecting ^' on 
these eigenvectors, 

<l^' = fii^»,<I>'W»„ ^^i' = (g^^\'^')g(J,1 fori = 1,2, 3 and J = 1,2. (3.19) 



The field ^^ is the neutral Higgs ghost field, while the fields ^^ and <I>^ form the 
charged Higgs ghost field, and the fields G^^' and G*-^-* together form charged Goldstone 
boson fields. 

If the vacuum expectation values broke the symmetry in even more directions, as in 
(3.1) with both U2 ^ and 5^0, the situation would be different: All 6 generators of 
SO{4:) would be broken, 4 of them corresponding to the 4 broken generators of the U{2) 
gauge group. Thus, there would still be 2 pseudo-Goldstone bosons. 

4. Concluding remarks 

In this paper we have analyzed the additional (approximate) symmetries which may arise 
in multi-Higgs-doublet models, due to the fact that the scalar potential may have more 
symmetries than required by the imposed gauge invariance. Moreover, for the kinetic terms 
(as a whole) we found that the symmetry group was U{N) x SU{k). In the case k = 2 
(i.e. scalar doublets) we found that the symmetry group of the kinetic terms, in the limit 
g' —7- 0, is enhanced to Sp{N) x 5*0(4). The most general C invariant Higgs potential (2.8) 

has the same SO (A) symmetry, only broken by the presence of the operator C^, that is, 

(3) (3) 

terms proportional to A^j, . In the case where A^^ is set to zero, we have a mass degeneration 

(3.12) (assuming vacuum alignment) between charged and C odd sectors in the limit g' — )■ 0. 

When we don't have vacuum alignment, but rather two broken (real) directions with the 

electromagnetic generator left unbroken, a pair of light, charged Higgs bosons should emerge 

(cf. section 3.3). 

A. P(A;,M), the symmetry group of (7^ 

We will here show that the set 

P{k, R) = {S£ GL2kimS^JS = ±J}, (A.l) 

given in eq. (2.17) is a Lie group: The associative law and the existence of the identity 
follow from GL2fc(M) (the set of all invertible, real 2k x 2k matrices) being a group. Define 

p-(fe, R) = {S€ GL2kiR)\S^JS = -J}. (A.2) 

The other component of P{k,M.) (what we could call P^{k,R)) is Sp{k,R). Then, if 
S~ G P~ and S~^ G Sp{k,R), then we easily see by the definition that 

S~S+,S+S- eP~{k,R), 

S+S+, S^S^ G Sp{k, R). (A.3) 
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So the set P(k,M.) is closed under group multiplication. This set also includes the inverse 
of each element. We only have to show this for elements S G X~ , since we already know 
Sp{k, M) is a Lie group. Let S^J'S = —J . Then 

{S^Y^S^JSS"^ = {S^)-^{-J)S-^, (A.4) 

and since we generally have that (^A^^^ = (^^^) , 

-J = {S-^fjS-^, (A.5) 

so 5"^^ G P^ too (still, P~ is not a group considered isolated, since it is not closed under 
group multiplication, and does not include the identity). 

We have now derived that P(A;, M) is a group. To prove it is a Lie group, we must prove 
that it is a (topologically) closed subset of GL2fc(lR): /(^) = A J' A is a continuous map, 
the set {±J} is closed in GL2fc(M), and hence P{k,R) = f~'^[{±J}] is closed in GL2fc(M). 

The determinant of P^{k, R) We will now show that the determinant of the matrices in 
the set P^(/c,M), consisting of the real matrices with the property S'^JS = —J, is (— 1)*^: 
First, we claim the set P~(/c,M) is given by 

p-{k,R) = Sp{k,R)C = eSp{k,R), (A.6) 

with C defined in eq. (2.13). This is so because if S' G P~{k,R.), then S'C G Sp{k,R) since 

{S'CYj{S'C) = C^{-J)C = J, (A.7) 

and then S' = SC for S* = S"C G Sp{k,M.), since C^ = /. Similarly with C S'p(A;,M). 
On the other hand, if S" G Sp{k,M.), then 

{SCfj{SC) = C^JC = -J, (A.8) 

so then SC G p-{k,R). Similarly, CS G P-{k,R). 

Now we can evaluate the determinant of an arbitrary element in S' G P~(/c,M). Since 
S' = SC for an element S G Sp{k, M), 

det(5') = det{S) det(C) = det(C), (A.9) 

since all matrices in Sp{k,M) have determinant 1 [20]. The determinant of a n x n matrix 
A can be written (sum over repeated indices) 

det(A) = e^i'--*"^i,,,---^„,,„ (A.IO) 

(the Leibniz formula). Then there is only one non-zero term in this sum for the matrix C, 
so the determinant is given by (no sum over k) 

det(C) = e''''-'^'Ci,iC2,2---C2k,2k = l'(-l)' = (-!)'• (A.ll) 

Hence by eqs. (A.9) and (A.ll), the matrices of P~{k,'K) have determinant (— l)'^. 
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Sp{k, M) and P^ {k, R) are not connected We want to show that Sp{k, M) and P^ {k, M) 
are two components of P{k,M.), i.e. they are not connected. Connected means the same 
as path connected for Lie groups. Assume that the two components are connected. Then 
there has to be a continuous path between e.g. / G Sp{k,M) and R G P^{k,R). Let X{t) 
be such a path, i.e. ^(0) = / and ^(1) = R, where X{t) is continuous. Consider the 
supremum 

to = svip{t\X^{t)JX{t) = +J}. (A.12) 

We know that X {!)"'" J' X (1) = —J . Moreover, consider the function 

f(t) = det{X^{t)JX{t) + J), (A.13) 

which is continuous for continuous functions X{t), since the determinant, matrix addition, 
muhiphcation and transposition are continuous. But f{t) is discontinuous for t = to, since 
there in any open interval contining to will be values t where f{t) = and other values 
where f(t) = det(2i7) = 2^^^, per definition of to- Hence our assumption that X(t) is 
continuous must be wrong, and hence the sets Sp{k,M.) and P~{k,M.) are not connected. 

B. The map p 

We will now introduce a map p which lets us easily translate between real and complex 
formulations of the kinetic terms we are studying. The map p preserves both matrix mul- 
tiplication, addition and the identity. ^^ We define p as a function from Mfc(C), the set of 
all A; X A; complex matrices, to M2fc(M), the set of all k x k complex matrices by 

With U a Lie group, p is a. Lie group isomorphism from U C Mfc(C) to p[U]. 

Now we want to show that the definition of p can be extended to vectors so that it 
preserve products of complex vectors and matrices: Let v be a complex k x 1 vector, let 
V = VFI + ivj, with vr, vj real and define 

^<"' ^ (':S) - (:) • <-' 

and 

p{v^) ^ (Re(^t), _Im(i;t)) = (v^, vj) . (B.3) 

Moreover, let A be a complex k x k matrix and let A = (Aji + iAj), with Ar, Aj real, 
then 

p{Av) = piA)piv), (B.4) 



^*p is an injective ring homomorphism [20]. On the other hand, the inclusion p[U{2)] C 50(4) shows 
that p does not preserve the determinant, even though it is a ring (or group) isomorphism on its image. 
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(Av)n\ _ (An -AA (vn 



since p{Av) = ( ^^^^^j = ( ^^ a i V") ~ Pi^)pi'^)- Furthermore, let u,v be complex kxl 
vectors, then 

Re{u^Av)= piu^)p{A)p{v), (B.5) 

since Re(t/t^t,) = Re[{u'^^-iuJ){AR+iAj){vR+ivi)] = («S uj) f^« '//) f^^) = p{u^ ) p{A) p{v) 
Then, 



Ai Ab. J \vi 



u^Av + v^A^u = p{u^)p{A)p{v) + p{v^)p{A^)p{u), (B.6) 

since the left hand side of eq. (B.6) equals its real part. 
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1. Introduction 

Future discoveries at the LHC may reveal a sector of scalar particles that is much richer than 
that of the Standard Model (SM). Some of the scalars may be responsible for generating 
the masses of fermions and the electroweak bosons [1], whereas others may be responsible 
for the dark matter [2, 3, 4, 5, 6, 7, 8]. Another interesting aspect of these theories is that 
CP violation is naturally accommodated, including its spontaneous breaking [9]. 

It is natural to classify such scalars according to their properties under the SU{2) 
associated with the electroweak sector of the Standard Model. In order to be compatible 
with electroweak precision data, one usually considers only SU{2) doublets and singlets. 
Even these representations are severely constrained by the data [10, 11, 12, 13]. 

Going beyond one or two doublets [14, 15, 16], one immediately has to face models 
having a large number of parameters. The structure of such potentials has been studied in 
[17]. 

Different doublets could be distinguished via their couplings to fermion fields. This 
idea is exploited in the so-called Model II version of the two-Higgs-doublet model (2HDM) 
[14], where one doublet couples to up- type fermions, and the other couples to down- type 
fermions. Another version of this idea is the one considered in Ref. [7], where each fermion 
or each family has its own Higgs field. 



We shall here consider instead the case when the different doublets can not be dis- 
tinguished (since we are not considering couplings to the fermions). An introduction of 
Yukawa couplings would naturally have broken the symmetry among the different dou- 
blets [18]. Thus, we shall here study the symmetries of models with N doublets — it turns 
out that by setting a certain S'0(4)-violating parameter A'^' to zero and assuming vac- 
uum alignment, the spectrum simplifies considerably. In particular, a certain "custodial" 
50(3) symmetry [19, 20] leads to a degeneracy between the mass matrix of the C odd (or 
equivalently, CP odd) and the charged Higgs bosons. 

This possibility of a symmetry group of the scalar polynomial which is larger than 
required by gauge invariance was pointed out by Weinberg [21] many years ago. In the 
theories considered the extra symmetry was assumed to be a symmetry of any quartic (i.e. 
renormalizable) potential of the scalar sector. 

The Standard Model with its single Higgs doublet is an example of a theory were the 
most general scalar potential has an extra 0(4) symmetry not contained in the SU(2) x U{1) 
gauge symmetry. An extension of the SM with an extra Higgs doublet (as for instance 
required if we want to introduce supersymmetry) or more, and with C symmetry, also has a 
quadratic potential which is automatically 0(4) symmetric. This is sufficient to enforce 0(4) 
mass relations up to perturbative corrections in the parameter g' . (Hence these corrections 
involve the gauge bosons Z and 7.) For the quartic potential, the extra 0(4) symmetry is 
broken by the parameter(s) A*-^-*, cf. eq. (2.8). Standard renormalizability instructs us to 
include into the Lagrangian, all terms allowed by the SU{2) x U{1) symmetry, hence it may 
be inconsistent to leave out terms proportional to X^^'. The presence of such parameters 
will in general lead to an order A*-^^ tree-level breaking of the additional symmetry. This is 
in principle not different from having the symmetry broken by loop corrections. It becomes 
in any case a question of the magnitude of the perturbation. 

We present a detailed classification of the possible terms in the potential, discussing 
charge conjugation and how the familiar custodial 0(4) symmetry of the SM potential gen- 
eralizes to an 0(4) (or, more generally, 0{2k)) symmetry for certain terms of the potential 
in the NHDM. The maximal 0{2k) symmetric potential turns out to be a more constrained 
potential than the maximal charge-invariant one. The kinetic terms are in general sub- 
ject to an independent classification, depending on the U{1) hypercharge coupling g' . We 
identify, in certain situations, a charged pair whose mass vanishes with g' . 

The paper is organized as follows. In section 2 we discuss the general potential and 
classify the corresponding Lagrangian, including the kinetic terms. In section 3 we dis- 
cuss spontaneous symmetry breakdown, and in section 4 we conclude. A couple of more 
matematical discussions are delegated to appendices. 



2. The NHDM potential and Lagrangian 

We define the iV-Higgs-doublet model, abbreviated NHDM, to be a system of A^ two-com- 
ponent complex scalar fields ^1,^2, ■■■ , ^N, each with the same transformation property 
under SU{2)l x U{1)y as the Higgs field of the Standard Model, and with dynamics defined 



by the Lagrangian density 

C{x) = Y, [D^'^m{x)]HD^<^^{x)] - y($i, $2, . . . , ^jv), (2.1) 

m 

where I^(<&i, <^2; • • • ) ^n) is a potential that — in its most general form — is given by (2.5) 
below. The covariant derivative D^ is defined as 



D^ = a^ + igTiW!" + ig'YBf", 



(2.2) 



where Wj^ and B^ are the SU{2)l and U{1)y gauge fields, respectively, and Tj = -^^i 
are the generators of SU{2) with 0"i,0"2,cj3 the Pauli matrices. Thus, our Higgs fields are 
labeled by two indices: The row index running from 1 to A^ is often written out explicitly as 
above, and an often hidden group index acted on by the gauge group. The latter are acted 
on by the matrices Tj in (2.2) (whose indices are also hidden). When written explicitly we 
shall use Greek letters from the beginning of the alphabet. 

To write the most general gauge-invariant potential in a renormalizable NHDM in a 
compact way, we introduce a set of linearly independent^ hermitian operators invariant 
under local SU{2)l x U{1)y transformations (this is a generalization of the approach for 
the 2HDM in [22]): 

Am = '^In'^m, 
Bran = Wm^n + ^i^m) = ^<^L^n) = Ba, (2.3) 



2 



t 'T. _ $t $„) = lui{<Pl<Pn) = Ca. 



Due to (anti-)symmetry under interchange of m and n we may impose the restriction that 
1 < m < n < N, and introduce indices a, b, . . . labeling such pairs. An explicit invertible 
encoding is 

1 < o = a{m, n) = m + -(n - 2)(n - 1) < -N{N -l)=Af. (2.4) 

We let m{a),n{a) denote the inverse of this encoding. We will use the summation convention 
that repeated indices from the beginning of the alphabet are summed from 1 to M, and 
repeated indices from the middle of the alphabet are summed from 1 to A^. The most 
general potential in the NHDM thus becomes a linear combination of all different quadratic 
and quartic factors in the $m (and <J>m) which can be formed from Am, Ba and C^:^ 



Fg($l, ...,<^n)= Aii'Mm + /"i'^^a + f^'i^Ca + A« 1„J„ + A^JS^S. 



b 



+ >^^SCaCb + X^J^lAmBa + X^lAmCa + Xffj BaCb, (2.5) 

where the "g" in Vg denotes "general". To avoid double counting we introduce the restriction 
m < n in the term involving Am4, and the restriction a < 6 in the terms involving X^b , 



^There are no linear relations between the operators in (2.3). However, they are algebraically dependent 
when N > 3, being restricted by (A'' — 2)^ polynomial equations of 8'th order in the fields. 

^Since (^la^^A ($J„o-^$„) = - (^l^i) i^L^n) + 2 (^l^„) {^l^^^i), other quartic invariants may 
be expressed by those chosen. 



AJjJj and A|j^ . We will not consider terms of degree higher than four, because these would 
destroy the renormalizability of the model [24]. From the hermiticity of the potential V^ 
all parameters // and A in the expansion (2.5) must be real. Thus the number of free real 
parameters in (2.5) is 

Ntot = N + 2J\f + ^N{N + I) + ATiAT + 1) + 2NAr + Af^ = ^N\N^ + 3), (2.6) 

which for A^ = 1 gives us the 2 parameters of the Standard Model (^u^ and A). N = 2 gives 
us the usual 14 parameters for the 2HDM. There are 54 parameters for N = 3 and 152 
parameters for A^ = 4. 

This counting does not take into account the fact that we may make SU{N) row 
transformations on the fields <J>m to eliminate some terms in (2.5). One possible choice 

(2) (3) 

is to transform the quadratic terms into a diagonal form, i.e. so that /ia = Ha =0. 
This in general leaves a matrix of A^ — 1 independent diagonal phase transformations (such 
that the determinant is unity). We may for instance use it to transform all A-[^ with 
n{a) = m{a) + 1 to zero. This reduces the number of parameters by N"^ — 1, i.e. to 
-^tot ~ 2 (-^^ ^ ^^ ^ ^)' yielding 11 parameters for N = 2 (in agreement with Davidson 
and Haber [15]), 46 parameters for A^ = 3, and 137 parameters for A^ = 4. 

2.1 The most general C-invariant NHDM-potential 

The charge conjugation operator C is a linear operator, multiplicative in the fields, which 
leaves complex constants unaltered, but maps fields onto their hermitian conjugate trans- 
posed; C{z^rn) = z^m , wherc z is a complex number.^ Then C{Ca) = —Ca, in con- 
trast to C{Am) = Am, and C{Ba) = Ba- We obtain a C-invariant potential by leaving 
out all terms which are odd in Ca, i.e., terms involving jia, ^ma, and A^^ . There are 
M + NJ\f + J\f'^ = \N{N - l)(Af2 + N + 2) such terms, leaving 

Nc = \n{N^ + bN + 2) (2.7) 

free parameters for the general renormalizable (7-invariant potential, 

Vb(^i, ...,^n) = li\k^Am + ^l^i^Ba + A« 1„A„ + AiJs.Sft 

+ \^^^CaC^ + \tiAmBa. (2.8) 

For A^ = 1 we get the usual 2 parameters of the standard model, the Higgs potential of which 
is automatically C-invariant. For A^ = 2 we get the usual (see e.g., [22]) 10 parameters. 
For A^ = 3 we get 33 parameters, and for A^ = 4 we get 86 parameters. 

This counting does not take into account that we may make 0{N) transformations on 
the row of fields $„ to eliminate some terms in (2.8). A natural choice is to transform 

(2) 

the quadratic terms to diagonal form, i.e. so that jXa =0. This reduces the number of 
parameters by \N{N — 1), i.e. to A"^ = \N{N^ + 3N + 4). This gives 9 parameters for 
N = 2 (in agreement with [16]), 30 parameters for A^ = 3, and 80 parameters for A^ = 4. 

The difference A^phases = A't'ot - ^c = j^^i^^ - 1) + 1 - A^ counts the number of 
genuine C-violating parameters in Vg (in agreement with Branco et al. [23]). 

^This definition assumes that we for some reasons have decided on a decomposition of all fields into their 
real and imaginary parts. It is not invariant under complex transformations of the fields, see e.g. [23]. 



2.2 Symmetries of A, B, C and C^ 

For generality we here consider k (rather than 2)-component fields, i.e. with SU{k) x U{1) 
as gauge group. To make it simpler to explore all possible symmetries we express the field 
^m in terms of its independent real (hermitian) components, <!>„ = ^m + i0m- Define 
2k X 2k matrices 

I=(^'^'], J=( ^'^'], (2.9) 

\0k Ik J ' \-Ik Ok J ^ ^ 

where the subscript k indicates the linear dimension of the submatrix involved. The complex 
scalar product between two fields $„ and <!>„, invariant under unitary {U{k)) transforma- 
tions, can be expressed in terms of two real bilinear forms^ 

Re(cl.t„<l>„) = B^n = (C 0^) ^ [ 0^ = <^n + 0™0n, (2.10a) 

iM'i'l^n) = C^n = (^^, 0™) ^ f Q J = C0n " ©l^n- (2.10b) 

The first is the Euclidean dot product between 2/c-component real vectors, the second is 
the Poisson bracket (symplectic form) of the same quantities viewed as coordinates and 
momenta of 2/c-dimensional phase space. The quantities in (2.10) are individually invariant 
under transformation groups larger than U{k). The first form B (with A as a special case) 
is invariant under the 0{2k) group of real orthogonal transformations, 

(:::)-o(*»). oro^i, ,2.n) 

the second form C is invariant under the Sp{k, R) group of real symplectic transformations, 

(::)-*(::). srjS^J. ,2.2) 

In this formulation the charge conjugation operator C discussed above can be represented 
as a particular 0{2k) transformation when acting on the fields ^„ and 0„ 

Considering infinitesimal transformations, O =X + eL + C(e^), S = X + eM + ©(e^), the 
conditions (2.11) and (2.12) become 

L^X + XL = 02fc, M^J + JM = Q2k. (2.14) 

Thus L must be a 2A; x 2k antisymmetric real matrix; there is a set (Lie algebra) so{2k) of 
2k'^ — k linearly independent such matrices. Writing out the condition for M in terms of 



'(^m,*™) = $^ • $„ = Bmn + iC„ 



k X k submatrices we find that it must liave tlie form 

^=(c-A^)' B = B^,C = C^. (2.15) 

There is a set sp(k) of A;^ + k{k + 1) = 2/c^ + k Unearly independent such matrices. The 
infinitesimal transformations of tlie original U{k) are the intersection of the sets so{2k) and 
sp{k). I.e., the matrices of the form ( r. 4 )) with A = —A and B = B. There are 
2 A;(A; — 1) + 2^1^ + 1) = /c^ such linearly independent matrices. 

The symmetries of C^ The form C^ (or more precisely CmnCm'n') has a bigger sym- 
metry group than the form C. Still, we will see that such operators (forms) will violate the 
full 0(4) symmetry we can assign the rest of the Lagrangian. In analogy with (2.12), C^ 
symmetries are given by 

(::)-*(::). s^js.^j, ,2.6) 

which can be collected in a set 

P{k, R) = {Se GL2kmS^JS = ±J}, (2.17) 

which we in Appendix A show is a Lie group. 
The component 

P'{k, R) = {S€ GL2kmS^JS = -J}, (2.18) 

consists of matrices with determinant 

det(p-(A;,IR)) = (-l)^ (2.19) 

as shown in appendix A. The group P{k,M.) = Sp{k,R) U P~{k,M.) will have the same Lie 
algebra as Sp{k,R), since the new component P"{k,M.) is not connected with the identity. 
This is manifested by the equation corresponding to eq. (2.14), 

J + e{M^J + JM) = ±J (2.20) 

not having any solution for the case of a —J' on the right side, see appendix A for a proof. 

The most general 0(2A:)-symmetric potential We can conclude that the most general 
0(2/c)-invariant potential can be written 

Vb(2fc)(^l, . . . , ^7V) = ^^^^m + /ui^^^a + A^ 1^1„ 

+ X^^BaB, + X^J^lA^Ba, (2.21) 



since we have seen that operators not containing any factor C are 0(2/c)-invariant. We 

(3) 
obtain an 0(2A;)-invariant potential by leaving out terms proportional to A^^ from the C 



invariant potential Vc (2.8). The number of terms in Vo(2fc) is then No{2k) = ^C ~ -^^ = 
\N{N^ + 5iV + 2) - \N'^{N - 1)2, giving 

Nom = \N{N + lf (2.22) 

free parameters for the general renormalizable 0(2/c)-invariant potential Vo[2k)- Fo^" -^ = 1 
we get the usual 2 parameters of the standard model, the Higgs potential being automati- 
cally 0(2A;)-invariant. For iV = 2 we get the usual 9 parameters, one parameter less than 
for the (7-invariant potential. For iV = 3 we get 24 parameters, and for iV = 4 we get 50 
parameters. 

This counting does again not take into account that we may make 0{N) transformations 
on the row of fields <!>„ to eliminate some terms in (2.21). We may transform the quadratic 

(2) 

terms to diagonal form, so that jXa =0. This reduces the number of parameters by 
\N{N - 1), i.e. to A^o(2fc) = \N{N'^ + N + 2). This gives 8 parameters for A^ = 2, 21 
parameters for A^ = 3, and 44 parameters for A^ = 4. 

2.3 Symmetries of the NHDM potential 

Since the NHDM-potential Vg is constructed from the invariants (2.3) the symmetries of the 
latter are reflected in the symmetries of the former, but in a manner depending on details 
of the construction: 

1. If Vg depends only on the A^'s, i.e. if only the parameters ix„l and \mn are nonzero, 
then the symmetry group of Vg is at least^ ®m=i 0{2k), since we can make indepen- 
dent transformations on each <!>„. 

2. If Vg depends only on the vim's and the BaH, i.e. for a C-invariant theory (2.8), where 
in addition the parameters A^^ = 0, then the symmetry group of Vg is at least 0{2k). 

(2) (2) 

It may contain several such factors if some of the parameters fia and A^^ vanish. To 

(2) 

analyze this we partition the <I*rn's into sets: If a parameter Ha is nonzero, then the 
fields ^m(a) and ^n(a) belong to the same set, with m(a) and n{a) denoting that m 

(2) 

and n are contained in a. If a parameter A^^ is nonzero, then the fields ^m(a) and 
^n(a) belong to the same set, and the fields ^m(b) and ^n(b) belong to the same set. 
With this partitioning into a maximal number of sets we may make one independent 
0{2k) transformation for each set. 

^ (3) (3) 

3. If Vg depends only on the Ca's, i.e. with only the parameters fia and A^^ being 

nonzero, then the symmetry group of V^ is at least Sp{k, M). If we (in the same manner 
as above) can partition the fields into several sets, then we may make independent 
Sp{k, M) transformations on fields belonging to different sets. However, since the 
additional symmetries in this case fail to be symmetries of even the zero'th order 
kinetic term (2.30), their significance is uncertain. 



^It could possibly be larger, since there might be additional row symmetries transforming fields $„ with 
different m into each other; such symmetries would require special relations among the parameters /ifn and 



A 



(1) 



4. With all parameters arbitrary the symmetry group of Vg is just the original SU{k) x 
C/(l) gauge symmetry. 

In this work we will pay special attention to the second scenario, with k = 2. 

2.4 Symmetries of the kinetic terms 

We now turn to the (global) symmetries of the kinetic terms of the Lagrangian, 

N 

K = Y,W + G''')<I>„(x)]t[(a^ + G^)^n{x)], (2.23) 

n=l 

with 

G^' = igTiWl" + ig'TBi". (2.24) 

Let Ki denote the terms of the i'tli order in the gauge fields. 

Consider the transformation of the kinetic terms linear in the gauge fields, Ki, under 
the map p defined in appendix B. We can then write^ 

N 

n=l 

N 



= Y,P{d^{^n)l)p{G,)p{i^n)k)+p{{^ni)p{G'^^)p{d,{^n)k) 
n=l 

N 

= J^ a^^r^c^n + <^l{-rnd,'^n (2.25) 

n=l 

where the subscript k in {^n)k indicates this is the usual complex Higgs fc-plet (in the case 
k = 2 the usual complex Higgs doublet), while <!>„ is the 2k dimensional real vector 

*"=(e:)' '■'■'■" 

where {^n)k = ^n + ^©m and where we also use eqs. (B.2), (B.3) and (B.6). 

In eq. (2.25) we have applied the transformation p on the gauge terms G^ defined in 

eq. (2.24), 

p (G^") = r^, p (G'^'f ) = -r^, (2.27) 

where T then reads 

~[gWi^ + g'YB^ gW^ i' ^'''"^ 



^Disregarding so-called Schwinger terms — here terms proportional to i[9;j(^(a;), 0(a;)] for a scalar field 
—or, alternatively, reasoning classically. 



with Wj!^ = ^j W-^ T?, summed over the set of real symmetric generators T^ of SU{k), and 
Wj" = ^"^iW-^T^, summed over the set of imaginary antisymmetric generators T". For 
k = 2 the two sets are respectively {2^^) 2*^^} ^^'^ {2^^}' 

Finally, we consider the kinetic terms quadratic in the gauge fields, 



N 



K2 = Y.{<^n)lG^^G^{^n)k 
n=l 

N 

n=l 

N 

= -Y.'^nT^'^n. (2.29) 



ra=l 



The symmetries of Kq When we first ignore couplings to the gauge fields the remaining 
terms can be written 



N 



Ko = Y.12 (^''^"^ ^/^^«" + ^''®«" ^A^Q"") ' (2.30) 



n=l Q=l 



where the group index a labels the k components of $„ = ^„ + iB„. This term is invariant 
under rotation of all components {^na, ©no} into each other. I.e., the symmetry group of 
Kq is 0{2kN). The connected part of this group is SO{2kN), whose generators are all real 
antisymmetric matrices, Lmn,ai3 = ~Lnm,i3a (i.e. .^"^ = —L, where transposition refers to 
both sets of indices). 

The symmetries of Kq and Ki: Next, consider the terms linear in the gauge fields 
again, cf. eq. (2.25), 

^1 = E ((^.K^©^)^^ (q^) - {<^Ql)r'd, [lA\ (2-31) 

with T given in eq. (2.28). 

Consider now an infinitesimal transformation 5^m,a = Lmn,ai3^n,i3j ^ind T denoting 
the 2k X 2k antisymmetric matrix in equation (2.28) (in group indices a, /3 — in addition it 
is proportional to the N x N unit matrix in row indices). The requirement that this is an 
infinitesimal symmetry transformation for Ki is that L T+T L = 0. Or, when we restrict 
L to be antisymmetric so that it also is an infinitesimal symmetry transformation for Kq, 

J^mn,afi Ifi'y ~ lafi ^mn,l3"/ = 0. [2.32) 

In order to determine the allowed structure of L, we expand these matrices into terms of 
definite symmetries {L^^> symmetric, and L^'^' antisymmetric) in the mn indices: 



with S (symmetric) and A (antisymmetric) restricted by the constraint (2.32). The sum 
runs over all possible combinations of allowed matrices . We next note that the antisym- 
metric matrices T can be expanded in the set 



T = {Tti,T:j,j} = {r^^ Z], L: \A : "W- (2-34) 




By substituting (2.33) into (2.32) we are led to search for the set of 2k x 2k real matrices S 
and A which commute with T for arbitrary values of the fields W^^ and i?^. It is sufficient 
to verify that this property holds for all elements of the set T. Let 

^ = W ^^r\ ' ^ = -5' o^ ^ = ^- (2.35) 

Requiring commutativity [see eqs. (2.32) and (2.33)] with the three types of matrices in T 
we obtain the conditions 















(2.36a) 


iiij = Jj A22, A22JJ = 


■■T!x^^, 


X12T! = 


-Ti X21, 


X21T! = 


-T!X,2, 


(2.36b) 


11 = -^^22 5 -^12 = — -^^21- 












(2.36c) 



Using (2.36c) we find that Xn and X12 must commute with all matrices T^, Tf (assumed 
to form an irreducible representation). By Schur's lemma they must then be proportional 
to the kxk unit matrix, so that S (xX and A (x J . Thus, the Lie algebra of the symmetry 
group of Kq and Ki consists of elements of the form 

T, — T rM^I -A- 7 rM^I (2 371 

This is the Lie algebra of U{N) written in real variables. 

The symmetries of Kq and Ki in the limit 51' —t- 0: A more interesting situation 
arises if we remove J from the set T, as would apply to the limit g' — t- 0. Then we still find 
that Xii + X22 and X12 — -^1^21 must commute with all matrices T^, Tf, and hence must be 
proportional to the unit matrix. Further, Xn — X22 and X12 + X21 must commute with 
all Tf", but anticommute with all Tf. For k = 2, i.e. for the gauge group SU{2)l x U{1)y 
in the limit g' — )■ 0, we find that nonzero solutions of (2.32), 

Xii - X22 oc e = i a^, X12 + X21 oc e, (2.38) 

are possible [see eqs. (2.33) and (2.35)]. This means that the possible antisymmetric ma- 
trices A may be any linear combination of matrices from the set 

J \ , (2.39) 




■^Without the restriction (2.32) there would be \k{2k + l)N{N -l) + \k{2k-l)N{N + 1) = kN{2kN-l) 
independent terms, equal to the number of generators of SO{2kN). 
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where the 2x2 matrix e was defined in eq. (2.38). The set ^ is a basis of generators for 
SU{2). Thus, the Lie algebra of symmetry generators for Kq and Ki in this case consists 
of elements of the form 

Lmn=lL(^n+Y.^^^^n, (2-40) 

allowing all possible symmetric N x N matrices L^^' for each A. There are ^N{N — 1) + 
^N{N + 1) = 2N'^ + N independent terms, equal to the number of generators of the N x N 
quaternionic symplectic group Sp{N). The generators (2.40) generate Sp{N), where the 
elements of Q act as the quaternions i, j and k. 

The above results were again found under the assumption that the fields Wj^ are ar- 
bitrary, and kept constant under the transformation. Combined SU{k) transformations 
of the W^ and the <!>„ fields still remain a symmetry. This symmetry is enlarged to 
SU{k) X Sp{N) as g' -^ 0. (In the case g' ^ it is SU{k) x U{N).) In the case 
k = 2 and (7' — )■ the custodial 5*0(4) symmetry^ is contained in SU{2) x Sp{N) in 
the following way: 50(4) ^ 5^7(2)^ x SU{2)r C SU{2)l x Sp{N).^ (More presicely, 
50(4) ^ {SU{2)l X SU{2)r)/Z2 [31].) The group SU{2)r C Sp{N) is generated (through 
exponentiation) by choosing L^^> = QnxN and each L^^> oc InxN in eq. (2.40). The gener- 
ators of 50(4) are hence the 3 generators in Q plus the 3 generators of the SU{2)l gauge 
group (written in real form). Finally, the U{1)y hypercharge symmetry group is contained 
in SU{2)r [20]. 

In the more general case of an SU{k) x Sp{N) symmetry of Kq and Ki in the limit 
gf' — 7- there is, in the same manner as above, a "custodial" SU{k) x SU{2)r symmetry, 
which also will contain SU{2) x SU{2)r = 50(4) subgroups. 

The symmetries oi K2 Next, consider the terms quadratic in the gauge fields cf. eq. (2.29), 

N 



n=l 




^2 = -E«'0")^'U" • (2-41) 



As in the symmetry analysis of Ki we want to find all matrices X such that XT^ = T^X. 
All matrices X which commute with T will fulfill this criterion (since T^ can be expanded 
in a set which consists of products of all possible pairs of matrices from T). Therefore, the 
symmetries of Ki are also symmetries of K2 ■ 

3. Spontaneous symmetry breakdown 

In this section we return to the case of A; = 2, i.e. with SU{2)l x U{1)y as the gauge group 
and a row of A^ Higgs doublets <!>„• Note however that many of our considerations have 
straightforward generalizations to fc > 2. 

As for the Standard Model, the potential Vg of equation (2.5), or Vc of equation (2.8), 
may acquire its minimum at nonzero values of the scalar fields, (^)o = ^^^', where <I> 



*The custodial SO(4) symmetry cannot be extended to an 0(4) symmetry [29]. 

'We are grateful to H. Haber, J. P. Silva and P. Ferreira for pointing out a mistake at this point in a 
previous manuscript. 
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(without a lower index) refers to the whole set of fields ^^n- This point, ^^^' , will belong 
to one or more manifolds of equivalent minima related by the symmetries of the potential. 
One may use these symmetries to transform ^^^' to a particular form. A possible one is to 
require for ^[ ' that only its lowest real component is nonzero. This can always be achieved 
by an SU{2)l x f/(l)y gauge transformation. Next, the upper component of $2 can be 
made real by the remaining f/(l) gauge transformation which keeps ^\ unchanged. Then 
one has no gauge freedom left to change $n for n > 3. However, it was shown by Barroso 
et al. [25] that a sequence of unitary row transformations can shift the vacuum expectation 
values to the first two fields of the row only^*^ , for instance (when written in complex form) 

<)=f°V c^f=(^J,V <l>(0)=0forn>3, (3.1) 

with vi, U2, V2 and 5 real. The special case ti2 = 5 = is usually referred to as vacuum 
alignment, in which case we may also transform V2 to zero by an orthogonal row trans- 
formation involving only <&i and $2- This is known as the Higgs basis [26]. However, for 
other purposes it may be more convenient to adopt a "democratic" basis in which the lower 
component of all (or most) fields <&„ have a nonzero real expectation value. It is related 
to the Higgs basis by an orthogonal row transformation which preserves the form of C and 
C/(l) electromagnetic gauge transformations (the latter preserving the definition of electric 
charge) . 

Assume now the case of vacuum alignment and a potential Vq/a) which is 0(4) invari- 
ant. Then the existence of the vacuum expectation values $(°' will break the (explicit) 
symmetry down to 0(3), with the consequence that the Higgs boson particle spectra and 
other physical properties will organize themselves into multiplets of 0(3) (broken by per- 
turbative corrections in g'). The number of broken symmetry generators is 3 whether we 
consider the symmetry broken from U{2) to f/(l) or from 0(4) to 0(3); this leads to the 
existence of 3 Higgs ghosts and no extra (pseudo-) Goldstone bosons^-^. 

3.1 Mass-squared matrices 

To make these statements slightly more explicit, as needed for calculation of the zero'th 
order (in g and g') particle masses, we expand the potential around $^°' to second order. 
There are no first order terms since we are expanding around a minimum. The matrix of 
second derivatives is the mass-squared matrix M'^^^o . It is restricted by symmetries to have 



^''One may collect the quantities <l?mL (a = 1,2) into two A'^-component complex vectors i'^' and $"'. 

By a U{N) row transformation one may first rotate "I'm so that only the component <l>j^ is nonzero, with 

"i>j real. There is a group of U{N — 1) transformations preserving this condition; this may be used to 

- (2) ~ (2) ~ (2) "" (2) 

transform "l>m so that only the components $1 , $2 ^.re nonzero, with $2 real. One cannot do better 

due to the existence of four real U{N) invariant parameters in ll^'^'ll, Hi^'H, and ^^^'^ ^^^' . But there 

remains a U{N — 2) group of transformations preserving this condition which can be used for other purposes. 

For an SU{k) x U{1) gauge group one may generalize this procedure to k vectors ^^^\ j — 1 . . . k. 

^^This remains true for general values of fe > 2; a set of aligned vacuum expectation values will break 

U{k) to U{k — 1) and 0(2k) to 0(2fc — 1). The number of broken generators is 2fc — 1 in both cases. The 

situation is different if we have two broken real directions, as in equation (3.1) with U2 = but 5^0. Cf. 

section 3.3. 
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a block diagonal form in the group indices a, /?. We use coordinates where <I>^ 
is expressed in terms of four real fields, 



$r 



^m + i0r 



(Pml + Wm.2 
Vm + 'nm + i0m3 



$(0) 







*m + i0r 



(3.2) 



It is now convenient to represent these on real form as 



$. 



Or 



<Pm2 



(3.3) 



We have the expansion 



,(0) 



y($W + A$) = {V)o + 



d^v 



2 \d^mpd<^r, 



A$„,pA$„, + 0(A$3 



(3.4) 



where ^mp denotes one of the four possibilities (f>mi, ilm, 4'm2, 4'm3, and the subscript 
indicates that a quantity is evaluated at <5 = $(^'. Now a set of generators for S'0(4)^^ is 



•h 



J2 



J3 



/O o\ 

0-10 

1 
\0 0/ 

/O -l\ 




\l 0/ 

where Ji, J2, J3 will transform the vacuum expectation value ^^^\ while J4, J5, Jg leave it 
unchanged, cf. Eq. (3.3). In terms of these, the broken generators of the SU{2) x U{\)y 
gauge group, written in real form by the transformation p defined in (B.l), are 



Ja 



( 1 o\ 
-10 


\ 0000/ 
/o o\ 



1 

\0 -1 0/ 



J5 



Jfi 



/O o\ 
0-1 


\o 1 oy 

/O -1 o\ 

10 

\0 0/ 



(3.5) 



1 



1 



V^-(J2 + J5), ^^'^-(Jl + J4), ^(1 

and the unbroken U{1) (electromagnetic gauge) generator is 



a 



J3, 



-(l + a3)^J6. 



(3.6) 



(3.7) 



If now V is invariant under 

-A<f>mi, A?7„ 



'ml 



Ar]jn, A<f>. 



'm2 



-A<t>. 



'm2, 



A4>. 



ms 



A<P: 



'm,3) 



there can be no terms in (3.4) mixing the sets {A(j)i,A(p2} and {Ar]m,A(j)m3}- If V in 
addition is invariant under 



A(t>, 



ml 



A0m2, Ar]m -^ Ar]„ 



A<j). 



'm2 



-A(f)ml, A<p. 



'm3 



A<j). 



'mSj 



^Equivalently 0(4): 50(4) and 0(4) have the same Lie algebra and hence the same generators. 
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there can be no terms in (3.4) mixing Ac/j^i and /S.(j)m2-, and we must have 

aV \ / d'^V 



ml ^(Pnl / \ <^rm2 <J(Pn2 / q 



^ch,mn- (3-8) 



We refer to this as the charged mass-squared matrix. The transformations considered 
generate a Z4 subgroup of the f/(l) gauge group generated by Jg, assumed to be a symmetry 
of V . We have formulated it this way as a reminder that invariance under discrete subgroups 
may be sufficient to impose useful restrictions on the mass matrices. 
If V is invariant under C transformations, 

Ar?m -^ Arym, A(/)m3 -> -^4>m.2,-, (3.9) 

(irrespective of how we define C to operate on the charged sector, e.g. A(j)mi — > A(j)mi, 
A.(prn2 -^ — A(/>m2) there can be no terms in (3.4) mixing Ar]m. and AcpmS- Thus the neutral 
mass-squared matrix decomposes into two more blocks, a C even and a C odd one. 

If V in addition is invariant under the transformations 

Acpml -^ A(j)ml, Arjm -> Ar/„, A(/)„2 -> -A(/)m,3, A0„3 -> A(/)„2, (3.11) 

which generate a Z4 subgroup of the SO (2) symmetry group generated by J4, we obtain 
the relation 

^C-,mn = ^ch,mn- (3.12) 

This explicitly demonstrates mass degeneracy between the charged and the C odd sec- 
tors [27]. Especially, if the potential is 0(4)-invariant (2.21), that is, we have a C-invariant 
theory where also the parameters X^f^ = [the latter implies (3.11)]^^, the above symmetry 
criteria for mass degeneracy are valid. Moreover, since the renormalization is not changed 
when the Higgs fields acquire a vacuum expectation value [28], we do not get any mass 

renormalization counterterms from the quartic operators. So even though 0(4)-violating 

(3) 
quartic terms proportional to A^^ cannot alone be prohibited by any discrete symmetry 

imposed on the NHDM Lagrangian [29], they will not show up as counterterms when renor- 

malizing the masses. Hence, the mass degeneration (3.12) will only be broken by loop 

corrections involving gauge bosons, since we get an exact 5*0(3) symmetry when g' = 0. 

With g' ^ and hence with an approximate 50(3) symmetry, the mass diff"erences of the 

charged and OP-odd sectors will be of order 0{g"^) oc 0{e^). 

On the other hand, the 50(3) symmetry between CP-odd and charged sectors could 

also be broken by counterterms of the type A^^^ C^, even though these terms are set to zero 

in the original potential, if we are considering scattering processes and not mass relations. 



'For supersymmetric theories we typically have \^^ -^ 0. 
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3.2 The Higgs ghosts 

Let A$ be chosen so that $(*^^ +eA$ + C'(e^) is a family of minima related by the symmetry 
of the potential V, 

^ -y($(0)+eA$) = 0, (3.13) 



to first order in e. By differentiating this relation with respect to e and then setting e = 
we find 

which reflects the fact that the matrix M^^^o has zero eigenvalues with corresponding 
eigenvectors A<&„^. We may take the latter to be A*^^*' oc Ji^^^' for i = 1,2, 3. Normalized, 

A$« = iv„„ 0, 0, 0)^ /a, A<5(^) = (0, 0, Vm, Of /a, A$(^) = (0, 0, 0, v^f /a, (3.15) 

with a^ = Ylm''^rn- The massless excitations in these directions correspond to a triplet of 
Higgs ghosts. There will be A^ — 1 additional 5*0(3) triplets of excitations in directions 
orthogonal to the ghosts. They correspond to physical particles. There will also be A^ 
5*0(3) singlets, transforming evenly under C, corresponding to physical particles. In the 
case of A^ = 2, the triplet is {H^, H^ , A), whereas the singlets are h and H [14]. 

3.3 Non-aligned vacuum expectation values 

We have assumed vacuum alignment in much of the previous discussions of this section. 
The phenomenologically most realistic deviation from this case is that we have a situation 
with two (real) broken directions, as in (3.1) with ii2 = but 5 7^ 0. This corresponds 
to a situation which preserves the U{1) electromagnetic gauge symmetry, generated by 
Jg. Its corresponding definition of electric charge is preserved, but the C symmetry (or, 
here equivalently-'^'^ , CP symmetry) is spontaneously broken. In this situation, assuming 



^*For the Lagrangian (2.1), spontaneous C and CP violation (SCV and SCPV) are equivalent: By 
definition [30], CP (C) is broken spontaneously if (1) There is a transformation that can be physically 
interpreted as CP (C) and which keeps the Lagrangian invariant and (2) There is no transformation that 
can be physically interpreted as CP (C) which keeps both the Lagrangian and the vacuum invariant. 
(SCV => SCPV): Assume C is spontaneously broken, and implement P by the spatial reflection 

P$„(t,r)pt ^$„(t,--r). (3.16) 

Hence condition (1) of SCPV is satisfied. Next, assume condition (2) for SCPV is not satisfied. Then there 
is a CP transformation 

(CP)$„(t, r)(CT)^ = Ug^^i^it, -r), (3.17) 

which leaves both the Lagrangian and the vacuum unaltered. But then the C transformation given by 

C$„(i,r)Ct = [/^„$r(t,r), 

with (7m„ — C/mn will infer that condition (2) of SCV does not hold, since the spatial reflection (3.16) does 
not change the vacuum nor the physics of the Lagrangian. This is a contradiction, and hence also condition 
(2) of SCPV must hold. 
(SCPV => SCV): Conversely, if CP is spontaneously broken and condition (1) hence is satisfied by 

(CP)<i>„(t,r)(CP)t = [/^r$r(i,-r-). 
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the potential only has the SU{2) x ^(1) gauge symmetry, the three SU{2) generators are 
spontaneously broken, and we hence only get three Higgs ghosts and no (pseudo-) Goldstone 
bosons. 

When we only have an SU{2) x U{1) symmetry in the potential, operators of the type 
C are present in the quadratic or in the quartic part of the potential. Two of the excitations 
of the charged mass squared matrix are still massless, corresponding to the usual charged 
Higgs ghosts. ^^ 

In the case of non-aligned VEVs (or operators C present in the quadratic part of the 
potential) the C even and odd excitations generally mix to give a 2N x 2N mass squared 
matrix for the neutral particles. Here one of the excitations will be massless, corresponding 
to a neutral Higgs ghost. 

On the other hand, assume that V is invariant under 0(4) transformations. Then there 
are no operators of the type C present in the potential. The explicit 0(4) symmetry is now 
broken down to 0(2) ~ U{1), so that 5 generators are broken. As before, 3 of these will 
generate excitations which correspond to the Higgs ghosts; the remaining 2 will correspond 
to nearly massless charged pseudo-Goldstone^^ bosons (massless to zero'th order in g'). 

To analyze the situation we again write ^ = ^^^' + $' in terms of real fields, 

^m = ^^^ + (0ml,??m,(Am2,Xm) with $(^) = (0, -Um,, 0, Ifm) . 

J4 and J5 are now also broken by the vacuum expectation values. Acting with the broken 
generators on ^^^> one finds five eigenvectors of the mass matrix with zero eigenvalues. 



we can re-implement the matrix U as U above, and re-implement P as the trivial transformation (3.16), 
and hence condition (1) of SCV is satisfied. Assume condition (2) of SCV does not hold. Then we, in a 
similar manner as for the case (SCV=>SCPV), can let U'^^ — U'^ and define P as in (3.16), and hence 
SCPV does not hold either, which is a contradiction. Hence condition (2) of SCV is satisfied. 

'^^The relation between the complex and real formulations of the charged mass squared matrix is given 
by the map p of appendix B: Let M^ denote the N x N complex mass squared matrix, and let M^ be the 
corresponding 2N x 2N real matrix. They are related by 



as can be seen from eq. (B.6). Then p{M^) = M^. The matrix M^ is Hermitian, and hence has only 
real eigenvalues. It then follows from the definition of p that A is an eigenvalue of M^ if and only if A 
is an eigenvalue of M^. Moreover, if v is an eigenvector of M^, (Re«,Imi')^ and (Imii, — Reu)''" will be 
eigenvectors of M^ with the same eigenvalue A. Hence oppositely charged particles will have the same mass, 
although terms (fiim and <j)2n mix by the presence of operators Cmn in the quadratic part of the potential or 
by complex (i.e. non-aligned) VEVs, and hence violate C. (This mass degeneration is a consequence of the 
remaining C/(l) — 50(2) hypercharge symmetry, i.e. that the generator Je is unbroken.) The identification 
(3.8) assumed that ijf>im and i/>2n did not mix, i.e. that Im(Mc ) — 0, cf. (B.l). 

^^Pseudo-Goldstone bosons stem from broken generators of the extra 0(4) symmetry of the potential, 
while Higgs ghosts by definition are generated by the broken generators of the gauge symmetry (which of 
course is a symmetry of the whole Lagrangian). The pseudo-Goldstone bosons acquire small masses from 
radiative corrections, and are hence not massless to all orders of perturbation theory, like true Goldstone 
bosons. True Goldstone bosons are, in contrast, generated by the spontaneous breaking of a symmetry of 
a total Lagrangian, not only a potential. 
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A<I>W oc Ji^'^'^\ After normalization 

A$« = (t;„,0,0,0)^/a, 

A$(^) = (0,0,-t-„,0)^/a, 

A$(^) = (0, w^, 0, -?;„)^ /\/a2 + 62, (3.18) 

A$W = (o,o,w;„, 0)^/6, 

A$(^) = («;„, 0,0, 0)^/6, 

where a^ = ^^ w^ and 6^ = ^^ it;^. We see that ah eigenvectors except A^m are in the 
charged sector, hence the two pseudo-Goldstone bosons (i.e. hght Higgs bosons) mentioned 
above will be charged, as claimed. The eigenvectors in (3.18) are normalized, but they are 
not necessarily orthogonal to each other. Their nonvanishing inner products are 

m 

Here | sini?| > 0, since the vacuum expectation values by assumption are non-aligned. Thus, 
the orthonormalized eigenvectors corresponding to the Higgs ghosts can be written 

f(?) = ^=i= (u;„,0,-z;„,0)^ = ^^ A<^g^ + ^^== A<pg\ (3.19) 

V a + V a + V a + 0^ 

H^^^ = J_^ {0,Wm,0,-v„,f = A<!>g\ 

where Hm oc Gj<l>'''\ Gj denoting the SU{2) generators as given by the map (3.6). The 
two eigenvectors corresponding to the Goldstone modes are orthogonal to those above, 

G(^) = 1 \a ("a^^^) + cos T? A^^^)") - b f A$(^) - cos ^ A$(^) 

sin 'd V a^ + b'^ 






G(2) = ^ \-a (a^^^^ - cos 1? A^^^A + b (a^^^'> + cos 1? A$(^) 

sin T? V a^ + &2 L V / V 



(3.20) 



They have been orthonormalized. We note that the normalization constant becomes infinite 
in the limit of aligned vacuum expectation values, siwd — t- 0. We recall that the set 
|//(-'^\ i?^^)^ //(3)^ (^(1)^ (j(2)| g^j.g j^g|. numerical eigenvectors of the mass-squared matrix. 
The corresponding zero mode fields are the quantum fields obtained by projecting ^' on 
these eigenvectors, 

^g'^ = fi^W, <^') Hj^l <^^;^ = (g^^\ <^') G^fH for . = 1, 2, 3 and j = 1, 2. (3.21) 



The field <l> is the neutral Higgs ghost field, while the fields ^ and <I> form the 
charged Higgs ghost field, and the fields G^"^' and G^"^' together form charged Goldstone 
boson fields. 

If the vacuum expectation values broke the symmetry in even more directions, as in 
(3.1) with both ti2 7^ and 5 7^ 0, the situation would be different: All 6 generators of 
50(4) would be broken, 4 of them corresponding to the 4 broken generators of the U{2) 
gauge group. Thus, there would be 2 pseudo-Goldstone bosons also in this case. 
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4. Concluding remarks 

We have analyzed the additional (approximate) symmetries which may arise in multi-Higgs- 
doublet models, due to the fact that the scalar potential may have more symmetries than 
required by the imposed gauge invariance. Moreover, for the kinetic terms we found that 
the symmetry group is SU{k) x U{N). In the limit g' ^ the symmetry group of the 
kinetic terms is enhanced to SU{k) x Sp{N), which has an SU{k) x SU{2) subgroup. In 
the case k = 2 the latter is the SU{2)l x SU{2)ji = 50(4) custodial symmetry. The most 
general C invariant Higgs potential (2.8) has the same 5*0(4) symmetry, only broken by the 
presence of the operator C^, that is, terms proportional to AJ^^ . In the case where X^f^ is set 
to zero, we have an exact mass degeneration (3.12) (assuming vacuum alignment) between 
charged and C odd sectors in the limit 5' — )■ 0. When there is no vacuum alignment, but 
rather two broken (real) directions with the electromagnetic generator left unbroken, a pair 
of light, charged Higgs bosons emerge (cf. section 3.3). 

The introduction of Yukawa couplings could further constrain the theory. With N 
doublets, one could imagine "simplified" models analogous to Model I and Model II for the 
2HDM, where only one doublet couples to all fermions, or where some doublets couple to 
up-type quarks, with others coupling to down-type fermions. Furthermore, with three or 
more doublets, one could arrange to let each fermion generation couple to its own doublet. 

If n doublets are without any Yukawa couplings, for example due to a discrete Z2 
symmetry, 

<^i^-^i, i = l,...,n (4.1) 

then such a sector would provide a dark matter candidate [5]. Indeed, with n > 1, there 
would be a whole "family" of states in this "inert" sector, some of which would carry electric 
charge. Those would therefore be observable. 

A. P{k,M.), the symmetry group of (7^ 

We will here show that the set 

P{k,R) = {Se GL^kimiS^JS = ±J}, (A.l) 

given in eq. (2.17) is a Lie group: The associative law and the existence of the identity 
follow from GL2A;(M) (the set of all invertible, real 2k x 2k matrices) being a group. Define 

P'{k,R) = {Sg GL2kiR)\S'^JS = -J}. (A.2) 

The other component of P{k,M) (what we could call P~^(k,M)) is Sp{k,M). Then, if 
S~ G P~ and S^ G Sp{k,M), then we easily see by the definition that 

S-S+,S+S- eP-{k,R), 

S+S^, S^S^ e Sp{k, M). (A.3) 



18 



So the set P(k,M.) is closed under group multiplication. This set also includes the inverse 
of each element. We only have to show this for elements S G X~ , since we already know 
Sp{k, M) is a Lie group. Let S^J'S = —J . Then 

{S^Y^S^JSS"^ = {S^)-^{-J)S-^, (A.4) 

and since we generally have that (^A^^^ = (^^^) , 

-J = {S-^fjS-^, (A.5) 

so 5"^^ G P^ too (still, P~ is not a group considered isolated, since it is not closed under 
group multiplication, and does not include the identity). 

We have now derived that P(A;, M) is a group. To prove it is a Lie group, we must prove 
that it is a (topologically) closed subset of GL2fc(lR): /(^) = A J' A is a continuous map, 
the set {±J} is closed in GL2fc(M), and hence P{k,R) = f~'^[{±J}] is closed in GL2fc(M). 

The determinant of P^{k, R) We will now show that the determinant of the matrices in 
the set P^(/c,M), consisting of the real matrices with the property S'^JS = —J, is (— 1)*^: 
First, we claim the set P^(/c,M) is given by 

p-{k,R) = Sp{k,R)C = CSp{k,R), (A.6) 

with C defined in eq. (2.13). This is so because if S' G P^(A;,IR), then S'C G Sp{k,R.) since 

{S'Cfj{S'C) = C^[-J)C = J, (A.7) 

and then S' = SC for S = S'C ^ Sp{k,R), since C^ = I. Similarly with CSp{k,R). 
On the other hand, if S" G Sp{k,R), then 

(SCfjiSC) = C^JC = -J, (A.8) 

so then SC G P-{k,R). Similarly, CS G p-{k,R). 

Now we can evaluate the determinant of an arbitrary element in S' G P~{k,R). Since 
S' = SC for an element S G Sp{k, R), 

det{S') = det{S) det{C) = det(C), (A.9) 

since all matrices in Sp{k,R) have determinant 1 [32]. The determinant of an n x n matrix 
A can be written (sum over repeated indices) 

det(A) = e^i'--*"^i,,,---^„,,„ (A.IO) 

(the Leibniz formula). Then there is only one non-zero term in this sum for the matrix C, 
so the determinant is given by (no sum over k) 

det(C) = ei'2'-'2fcCi,i C2,2 • • • C2k,2k = l'(-l)' = (-1)'. (A.ll) 

Hence by eqs. (A.9) and (A.ll), the matrices of P^{k,R) have determinant (— l)'^. 
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Sp{k, M) and P^ {k, R) are not connected We want to show that Sp{k, M) and P^ {k, M) 
are two components of P{k,M.), i.e. they are not connected. Connected means the same 
as path connected for Lie groups. Assume that the two components are connected. Then 
there has to be a continuous path between e.g. / G Sp{k,M.) and C € P^{k,R). Let X{t) 
be such a path, i.e. ^(0) = / and X{1) = C, where X{t) is continuous. Consider the 
supremum 

to = svip{t\X^{t)JX{t) = +J}. (A.12) 

We know that X {!)"'" J' X (1) = —J . Moreover, consider the function 

f(t) = det{X^{t)JX{t) + J), (A.13) 

which is continuous for continuous functions X{t), since the determinant, matrix addition, 
muhiphcation and transposition are continuous. But f{t) is discontinuous for t = to, since 
there in any open interval containing to will be values t where f{t) = and other values 
where f(t) = det(2i7) = 2^^^, by definition of Iq. Hence our assumption that X(t) is 
continuous must be wrong, and hence the sets Sp{k,]B.) and P~{k,M.) are not connected. 

B. The map p 

We introduce a map p which lets us easily translate between real and complex formulations 
of the kinetic terms we are studying. The map p preserves both matrix multiplication, 
addition and the identity.^ We define p as a function from Mfe(C), the set of all k x k 
complex matrices, to M2A;(M), the set of all 2k x 2k real matrices by 

With U a Lie group, /? is a Lie group isomorphism from U C Mfc(C) to p[U]. 

Now we want to show that the definition of p can be extended to vectors so that it 
preserves products of complex vectors and matrices: Let v be a complex k x 1 vector, let 
V = vr + ivj, with vr, vj real and define 

^<"' ^ (':S) - (:) • <-' 

and 

p{v^) ^ (Re(^t), _Im(i;t)) = (v^, vj) . (B.3) 

Moreover, let A be a complex k x k matrix and let A = (Aji + iAj), with Ar, Aj real, 
then 

p{Av) = piA)piv), (B.4) 



^^p is an injective ring homoniorphism [32]. On the other hand, the inclusion p[U{2)] C 50(4) shows 
that p does not preserve the determinant, even though it is a ring (or group) isomorphism on its image. 
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(Av)n\ _ (An -AA (vn 



since p{Av) = ( ^^^^^j = ( ^^ a i V") ~ Pi^)pi'^)- Furthermore, let u,v be complex kxl 
vectors, then 

Re{u^Av)= p{u^)p{A)p{v), (B.5) 

since Re(t/t^t;) = Re[{u^^-iuJ){AR+iAj){vR+ivi)] = ^ uj) (f^ 'fj) f^^) = p{u^ ) p{A) p{v) 
Then, 



Ai Ab. J \vi 



u^Av + v^A^u = p{u^)p{A)p{v) + p{v^)p{A^)p{u), (B.6) 

since the left hand side of eq. (B.6) equals its real part. 
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